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ensemble 


i8 








L/a 


B35.32 


1.95 


0.0035 


0.135 


0.170 


32 


B55.32 


1.95 


0.0055 


0.135 


0.170 


32 


B75.32 


1.95 


0.0075 


0.135 


0.170 


32 


B85.24 


1.95 


0.0085 


0.135 


0.170 


24 


D20.48 


2.10 


0.0020 


0.120 


0.1385 


48 


D30.48 


2.10 


0.0030 


0.120 


0.1385 


48 



Table 1: The ensembles used in this investigation. The notation of ref. is used for labeling the ensembles. 

1. Introduction 

The decay constants of the Pion, Kaon, the D- and the Dy-meson are phenomenologically 
interesting quantities, not least because the ratio fK/fn together with the well known \Vud\ can be 
used to determine \Vus\- While experimentally /;;: and fx are well known and fo and /d^ less so, 
lattice QCD is in principle able to provide calculations of all of these from first principles. And 
recent advances in the field allow now also for statistically precise determinations with Nf = 2 and 
Nf = 2 + \ dynamical quark flavours ^ ||, ^ ^, with systematic uncertainties more or less under 
control. 

In this proceeding contribution we are going to present yet another determination of the afore- 
mentioned decay constants, however with a dynamical charm quark in place, i.e. for QCD with 
Nf = 2 + \ + \ quark flavours. 

2. Set-up 

We use gauge configurations as produced by the European Twisted Mass Collaboration (ETMC) 
with Nf = 2 + \ + \ flavours of Wilson twisted mass quarks and Iwasaki gauge action. The set-up 
is described in ref. and the ensembles used in this investigation are summarised in table H} The 
twisted mass Dirac operator in the light - i.e. up/down - sector reads [Q] 

Di=Dv/ + mo + iHiYs T3 (2. 1) 

and in the strange/charm sector [Q] 

D/, =Div + mo + /AtCT75Ti +/I5T3, (2.2) 

where Dw is the Wilson Dirac operator. The value of mo was tuned to its critical value as discussed 
in refs. ^ in order to realise automatic iff (a) improvement at maximal twist [p^. Note that the 
bare twisted masses iJ.(j § are related to the bare strange and charm quark masses via the relation 

m,,, = AtcT ± (Zp/Zs) Ms (2.3) 

with pseudoscalar and scalar renormalisation constants Zp and 2$. 

In the strange and charm quark sector, we use a mixed action approach with Osterwalder-Seiler 
(OS) valence quarks by formally introducing a twisted doublet for valence strange and charm quark 
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j3 = 1.95 jS =2.10 

'a^s 0.0130 0.0110 

0.0145 0.0120 
0.0160 0.0130 
0.0180 0.0150 
0.0210 0.0180 

0.0200 

Table 2: Bare values of the valence quark masses in the strange region for /3 = 1 .95 and j3 =2.10. 



flavours Ql 1| , |12| ], a set-up without flavour mixing artifacts in the valence sector. The two actions 
in the sea and the valence sector can be matched by tuning the bare values of valence strange and 
charm such that unitary Kaon and D-meson masses are reproduced. Unfortunately, it turned out 
that the unitary Kaon and D-meson masses were not exactly tuned to their physical values and 
currently we are still lacking ensembles which would allow to interpolate to the corresponding 
physical values. Hence, for the time being we vary the valence quark masses to interpolate the 
valence Kaon and D-meson masses to their physical values, as discussed below. This approach has 
been successfully applied to the ETMC Nf = 2 flavour gauge configurations in ref. [Q]. 



The determination of the (unitary) K-meson mass is described in ref. []13|]. Its decay constant 
can be determined from 

A = (m, + m.)^^^, (2.4) 
where Pk represents the projection to the physical Kaon interpolating operator as discussed in 



ref [13]. iris is the strange quark masses defined in eq. (2.3). The pseudoscalar decay constant /ps 



in the valence sector is determined from 

/ps-(Mv«/+M™/J sinhmps ' ^^-^^ 
where P = 75^2 with quark fields q\ , ^2 suitably chosen for the desired quark content. The meson 



mass mps and the matrix element KOlPj/'i')! entering eq. (2.5) have been extracted from a single 



state fit of the corresponding two-point pseudoscalar correlation function. The replacement of mps 



with sinh mps in the lattice definition ( |2.5[ ) of the decay constant helps in reducing discretisation 
errors for heavy meson masses [Q]. 

Any value for /ps and mps depends on the mass-value of the light dynamical quark (sea strange 
and charm quark mass values are fixed for a given j8-value) and two valence quark masses, which 
we denote by fps{lJ.f,IJ-i"' ,1^1^"^)- We have investigated several values for the valence quarks: 
jj.™' = = He and five to six values in the strange and charm quark region /i™' = /iv.c- The 
ensembles used are summarised in table |l| and the bare values of the valence strange and charm 
quark masses in tables I and |. 

3. Results 

One may expect large differences in between the unitary set-up and the mixed action set-up at 
finite values of the lattice spacing. In order to investigate this point we have determined /k in the 
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alXc 



i = 1.95 
0.200 
0.215 
0.240 
0.260 



j3 =2.10 
0.1650 
0.1800 
0.2000 
0.2250 
0.2500 



Table 3: Bare values of the valence quark masses in the charm region for j3 = 1 .95 and j3 =2.10. 



unitary and the mixed action set-up after matching the Kaon mass. The result is shown in figure ^ 
where we plot /k as a function of the squared pion mass. Both quantities are in units of /o, the pion 
decay constant in the chiral limit as determined in ref. [^. Within errors unitary and mixed set-up 
determination of fx agree. Moreover, results for the two available values of the lattice spacing 
agree within errors. This indicates small lattice and small unitarity breaking artifacts at least in 
fx- Note that the ratio Z^/Z^ used for the determination of the unitary fx has been determined by 



matching the bare value of the OS strange quark mass to pia and /ig via eq. ( |2.3| ). 

In order to determine the physical value of we fit SU(2) ;^PT formulae to our pion and 
Kaon decay constants data [0, ^ simultaneously according to 



/ps(At£,jU^At.v) 



/o-(l-2^„ln^/,+Z.^/,) 

(/d ' "I" fm ^ss) ■ 

3 



l-^^^ln^z + C^'^^ + C^U^H 



where 



(3.1) 



(3.2) 



(3.3) 



(47r/o)2 

are expressed in our analysis as a function of meson masses^ . We correct our data for finite size 
effects using NLO x^T [ 



16] 



fps{l^£,l^e,l^e;L)=fpsil^e,He,l^e) ■ [I -2^iigi{L,^ii)] 



fps (At£ ,He,l^s;L)= fps Ms- ) 



(3.4) 



We fit to our j8 = 1.95 and j8 =2.10 data simultaneously. We do not include lattice artifacts of 
^(c?^) in our fit, because the amount of data is not sufficient to do so: including these effects lets 
the fits become instable. Moreover, the data is fittable without these terms with ;^^/dof = 50/30. 

The physical input to our fits is = 130.7 MeV, m,^ = 135 MeV and mk = 497.7 MeV. In 
figure ^ we show pion and Kaon decay constants as a function of the squared pion mass The 
corresponding Kaon decay constant is determined at a value of mp^{}jL^,}JLs:i^s) = 2m^ — m\. For 
the figure we have interpolated our data linearly to these values. As a result we obtain fx/ fn = 
1.224(13), fx = 160(2) MeV and I4 = 4.78(2) with statistical errors only as determined from a 
bootstrap analysis. 



We use the normalization in which /^j = 130.7MeV. 
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Figure 1: fKl fo as a function of the pion mass squared. Results for the unitary and the mixed action fK_ are 
shown for two values of the lattice spacing. 



The fit also allows us to determine the values of the lattice spacings at j3 = 1.95 and j8 =2.10 
and the value of /q. All these quantities agree very well with the results obtained in ref. 

Following the procedure described in ref. [H] we have also analysed the data for /o and /^^ 



using SU(2) heavy meson chiral perturbation theory [|l7|]. We consider here the expansions of 



fo^Jm), and 



(3.5) 



including terms proportional a^m^^ and 1 /wid,- For details we refer to ref. [Q]. Our first results for 
fo and fo, are very encouraging, however, quoting results requires a better control of the system- 
atics involved in this investigation. 



4. Summary and Outlook 

We have presented the first determination of the pseudoscalar decay constants fx, fo and fo^ 
from lattice QCD with dynamical up, down, strange and charm quark flavours. We have used a 
mixed action approach with Osterwalder-Seiler valence quarks on a maximally twisted mass sea. 
The analysis indicates small lattice and unitarity breaking artifacts. The preliminary results are 
Jk//^ = 1.224(13) and fx = 160(2) MeV obtained with a SU(2) chiral perturbation theory fit. 



The errors are statistical only. Using the results of ref. [18] this translates to \Vus\ = 0.220(2) 
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Figure 2: Pion and Kaon decay constants as a function of m^. We show data for two values of the lattice 
spacing a — 0.079 fm and a = 0.060 fm, corresponding to j3 = 1 .95 and j3 =2.10. 

A comparison of our results for fx and fx / fn to results available in the literature shows that 
to the current level of accuracy there is no difference visible, neither to Nf = 2 flavour results [Q] 
nor to A^/ = 2 + 1 flavour results ||3|, ||, [T]]. 

A similar analysis has been performed for fa and fo,- In the case of these charmed quantities it 
will be in particular interesting to understand whether lattice artifacts proportional to a^m^ are small 
enough in our lattice set-up to allow for a precise determination of the corresponding quantities. 
However, the Nf = 2 results presented in ref. give rise to optimism that also fj) and /o^ can be 
reliably determined in our set-up. 

Clearly the results presented here need a better understanding of the systematic uncertainties. 
This includes in particular the dependence of the decay constants on the sea strange and charm 
quark mass values. ETMC is currently producing ensembles to investigate this point. These new 
ensembles should eventually allow to interpolate to the physical values of mx and mo- Moreover, 
more ensembles with different light quark mass values are required for j8 = 2.10 to better control 
lattice artifacts, for which also a third value of the lattice spacing is desirable. 

On the analysis site we are currently implementing different fit formulae (see for instance 
refs. [^, l^]) in order to better understand the extrapolation in the various quark masses. 
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